WE CLAIM: 

1 . A method for determining one or more signal parameters for a plurality of 
tones in an input signal, the method comprising: 

5 a) receiving samples of the input signal, wherein the input signal includes a 

plurality of tones; 

b) generating a frequency transform of the samples, wherein the frequency 
transform comprises frequency response data for the plurality of tones; 

c) determining one or more signal parameter estimates characterizing each of 
10 the plurality of tones; 

d) for each respective tone of the plurality of tones, generating modified 
frequency response data for the respective tone using one or more signal parameter 
estimates of one or more of the other tones, wherein the modified frequency response data 
for the respective tone includes reduced effects of other tones of the plurality of tones; 

15 wherein step d) is performed in an iterative fashion; and 

e) determining final one or more signal parameter estimates for each of the 
plurality of tones based on the modified frequency response data of the plurality of tones. 

2. The method of claim 1, wherein step c) comprises: 

20 cl) determining a plurality of potential local peaks based on the frequency 

response data, wherein each potential local peak corresponds to a tone; 

c2) determining one or more signal parameter estimates characterizing each of 
the pluraHty of tones based on said plurality of potential local peaks; 

c3) storing said one or more signal parameter estimates for each of the 
25 pluraUty of tones; 

3. The method of claim 2, wherein step d) comprises: 
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for each respective tone characterized by said one or more signal parameter 
estimates, 

dl) removing effects of other tones of the plurality of tones from the 
frequency response data of the potential local peak of the respective tone using said one 
5 or more signal parameter estimates of each of said other tones to generate modified 
frequency response data for the respective tone; 

d2) applying a single tone estimation method to the modified 
frequency response data, thereby generating a refined estimate of the one or more signal 
parameters of the respective tone; and 
10 d3) updating said one or more signal parameter estimates with the 

refined estimate of the one or more signal parameters of the respective tone; 

wherein, after said removing, said applying, and said updating for each tone, the 
one or more signal parameter estimates comprise refined estimates for each tone. 

15 4. The method of claim 3, wherein step e) comprises: 

performing steps dl) through d3) one or more times to generate final one or more 
signal parameter estimates for each of the plurality of tones. 

5. The method of claim 4, 
20 wherein said one or more signal parameter estimates characterizing the plurality 

of tones comprise a knowledge base of said signal parameter estimates; 

wherein said storing one or more signal parameter estimates comprises: 

storing a copy of said knowledge base, wherein said copy comprises a 
prior knowledge base; and 
25 wherein said updating said one or more signal parameter estimates with the 

refined estimate of the one or more signal parameters of the respective tone comprises: 

updating said knowledge base with said refined estimate of the one or 
more signal parameters of the tone, wherein, after said removing, said applying, and said 
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updating for each tone, the knowledge base is comprised of said refined estimates for 
each tone. 

6. The method of claim 5, wherein said performing steps dl) through d3) one 
5 or more times comprises performing steps dl) through d3) in an iterative manner imtil 

one or more of 1) the number of iterations exceeds a threshold number of iterations, and 
2) the difference between values of the one or more signal parameters of said knowledge 
base and values of the one or more signal parameters of said prior knowledge base is less 
than a threshold value. 

10 

7. The method of claim 4, wherein said one or more signal parameter 
estimates comprise estimates for one or more of frequency, amplitude, and phase. 

8. The method of claim 7, wherein said estimate of frequency comprises the 
1 5 frequency value of the determined potential local peak. 

9. The method of claim 7, wherein said estimate of amplitude is zero. 

10. The method of claim 7, wherein said estimate of amplitude is the value of 
20 the potential local peak. 

1 1 . The method of claim 7, wherein said estimate of phase is zero. 

12. The method of claim 4, wherein said frequency transform comprises a 
25 windowed Fast Fourier Transform (FFT). 

13. The method of claim 4, ftuther comprising: 

storing the final signal parameter estimates for each of the plurality of tones. 
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14. The method of claim 4, further comprising: 

outputting the final signal parameter estimates for each of the plurality of tones. 

5 15. The method of claim 4, wherein said storing said one or more signal 

parameter estimates for each of the pluraUty of tones comprises: 

sorting said signal parameter estimates for the plurality of tones by amplitude of 
said potential local peaks; and 

wherein said removing effects, said applying a single tone estimation, and said 
10 updating said one or more signal parameter estimates are performed for each of said 
plurality of tones in order of decreasing amplitude. 

16. The method of claim 4, wherein said removing effects of other tones fi*om 
the frequency response data of the potential local peak of the respective tone comprises: 

15 determining two or more sample values proximate to the potential local peak; and 

for each of the other tones, 

calculating an effect value for each of the two or more sample values using 
the estimated signal parameters for the tone; and 

subtracting the respective effect values firom each of the two or more 
20 sample values to generate a corresponding two or more modified sample values, wherein 
said two or more modified sample values comprise said modified fi-equency response 
data; 

wherein said modified fi-equency response data are usable in applying said single 
tone estimation method. 

25 

17. The method of claim 16, wherein said two or more sample values are 
consecutive sample values. 
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18. The method of claim 1 6, 

wherein said frequency transfomi comprises a windowed Fast Fourier Transform 
(FFT) with window function W(f)^ wherein / denotes frequency; 

wherein said estimated signal parameters for the tone comprise an amplitude ^4/, a 
5 frequency //, and a phase (f>u and wherein said calculating an effect value for each of the 
two or more sample values using the estimated signal parameters for the tone comprises 
calculating: 

dF{m) = A.e^'^Wif^ + /,) + Ae^^W{f^ - /,) , 

10 

for each sample w, wherein m denotes an index of the sample, andy», denotes the 
frequency bin of the sample; and 

wherein subtracting the respective effect value from each of the two or more 
sample values comprises subtracting a respective value of dF(m) from each of the two or 
1 5 more sample values. 

19. The method of claim 4, wherein said single tone estimation method 

comprises: 

identifying a frequency location proximate to an amplitude peak in the frequency 
20 transform, wherein the amplitude peak corresponds to the tone; 

selecting two or more frequency bins proximate to the frequency location in the 
frequency transform; and 

determining a tone frequency value that minimizes a difference between at least a 
first expression and a second expression; 
25 wherein the first expression comprises a sum of two or more numerator 

terms divided by a sum of two or more denominator terms, wherein the first expression 
includes a tone frequency variable; 
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wherein each numerator term and each denominator term corresponds to 
one of the frequency bins; and 

wherein the second expression comprises a sum of one or more of the 
numerator terms divided by a sum of one or more of the denominator terms, wherein the 
second expression includes the tone frequency variable; 

wherein the tone frequency value comprises a frequency of the tone. 

20. The method of claim 19, 

wherein the tone frequency variable represents a correct tone frequency value of 
the tone; and 

wherein the determined tone frequency value represents an approximation of the 
correct tone frequency value. 

21. The method of claim 20, 

wherein the first expression is approximately equivalent to the second expression 
when the correct tone frequency value is used for the tone frequency variable in the first 
and second expressions. 

22. The method of claim 19, 

wherein a ratio of each numerator term and its corresponding denominator term 
represent an amplitude of the tone at a respective bin. 

23. The method of claim 19, 

wherein a ratio of each numerator term and its corresponding denominator term 
represent a complex amplitude of the tone at a respective bin. 

24. The method of claim 19, 

wherein said determining a tone frequency value comprises: 
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computing a plurality of differences between the first expression and the 
second expression for different respective tone fi-equency values of the tone fi-equency 
variable; and 

selecting the tone frequency value that produces a smallest difference. 

5 

25. The method of claim 24, wherein said computing a plurality of differences 
comprises performing a Newton-Rhapson root finding method. 

26. The method of claim 19, 

10 wherein the first expression and the second expression are each a complex 

expression. 

27. The method of claim 26, 

wherein said determining a tone frequency value that minimizes a difference 
15 between a first complex expression and a second complex expression comprises 
determining a tone frequency value that minimizes a difference between an amplitude of 
the first complex expression and an amplitude of the second complex expression 

28. The method of claim 19, wherein said frequency transform comprises a 
20 windowed Fast Fourier Transform (FFT), and wherein the first expression and the second 

expression have the form: 



tXnn)W\f„ -f^)-F'(n)W(f„ +/.)) 

k+\ 



Z{nL-f>)f-\nL+f>)f) 



A+1 



I 

A+1 



25 



wherein: 

F(n) is the wth value of the single sided scaled FFT spectrum; and 
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W represents a window function, wherein the window function is shifted 
by a value of the tone frequency variable fj. 

29. The method of claim 19, wherein said frequency transform comprises a 
windowed Fast Fourier Transform (FFT), and wherein said generating a frequency 
transform of the samples comprises generating a power spectrum of the samples; 

wherein the first expression and the second expression have the form: 

k+M k+M-l 

j:i\F(n)\x\w(f„ -/)|) -/•)!) 

k+\ _ k+\ 

k+M k+M-\ 

B^(fn-fi)\ mnL-M 

k+\ *+i 

wherein: 

F(n) is the nth value of the single sided scaled FFT spectrum; and 
W represents a window function, wherein the window function is shifted 
by a value of the tone frequency variable f^. 

30. The method of claim 19, further comprising: 

computing one or more of the amplitude and phase of the tone using the 
determined tone frequency value. 

31. The method of claim 1 9, 

wherein at least one of the two or more frequency bins is on each side of the 
frequency location. 

32. The method of claim 1 9, 
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wherein said identifying comprises identifying a frequency location of at least one 
first magnitude peak in the frequency transform; and 

wherein said selecting comprises selecting two or more frequency bins proximate 
to the at least one first magnitude peak in the frequency transform; 

5 

33. The method of claim 4, wherein the frequency transform comprises a 
frequency transform array for the plurality of tones, and wherein said single tone 
estimation method comprises: 

identifying a frequency location of at least one first magnitude peak in the 
10 frequency transform array; 

selecting two or more frequency bins in a neighborhood of the at least one first 
magnitude peak in the frequency transform array; and 

determining a tone frequency value that minimizes a difference between a first 
expression and a second expression; 
15 wherein the first expression comprises a sum of two or more numerator 

terms divided by a sum of two or more denominator terms, wherein the first expression 
includes a tone frequency variable that represents a correct tone frequency value of the 
tone; 

wherein each numerator term and each denominator term corresponds to 
20 one of the frequency bins; and 

wherein the second expression comprises a sum of one or more of the 
numerator terms divided by a sum of one or more of the denominator terms, wherein the 
second expression includes a tone frequency variable that represents the correct tone 
frequency value of the tone; 
25 wherein the determined tone frequency value represents an approximation of the 

correct tone frequency value. 

34. The method of claim 33, 



Atty. Dkt. No.: 5150-62101 



Page 53 



Meyertons, Hood, Kivlin, Kowert & Goetzel PC 



wherein the first expression is approximately equivalent to the second expression 
when the correct tone frequency value is used for the tone frequency variable in the first 
and second expressions. 

5 35. The method of claim 33, 

wherein a ratio of each numerator term and its corresponding denominator term 
represent a complex ampHtude of the tone at a respective bin. 

36. The method of claim 4, wherein said single tone estimation method 

comprises: 

identifying two or more frequency bins proximate to a first magnitude peak in the 
frequency transform; and 

selecting a tone frequency value that makes a plurality of expressions most nearly 

equal; 

wherein each of the plurality of expressions comprises a sum of one or 
more numerator terms divided by a sum of one or more denominator terms, wherein each 
of the pluraUty of expressions includes a tone frequency variable; 

wherein each numerator term and each denominator term corresponds to 
one of the frequency bins; and 

wherein a ratio of each numerator term and its corresponding denominator 
term represent a complex amplitude of the tone at a respective bin; and 

wherein the tone frequency value comprises a frequency value of the tone. 

37. The method of claim 36, 

25 wherein the tone frequency variable represents a correct tone frequency value of 

the tone; and 

wherein the determined tone frequency value represents an approximation of the 
correct tone frequency value. 



10 



15 



20 
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38. The method of claim 37, 

wherein the first expression is approximately eqxiivalent to the second expression 
when the correct tone fi-equency value is used for the tone frequency variable in the first 
and second complex expressions. 



10 



39. The method of claim 36, 

wherein said determining a tone frequency value comprises: 

computing a plurality of differences between the first expression and the 
second expression for different respective tone frequency values of the tone frequency 
variable; and 

selecting the tone frequency value that produces a smallest difference. 



40. The method of claim 36, 

15 wherein a ratio of each numerator term and its corresponding denominator term 

represent a complex amplitude of the tone at a respective bin. 

41. The method of claim 36, wherein said frequency transform comprises a 
windowed Fast Fourier Transform (FFT), and wherein the first expression and the second 

20 expression have the form: 



k+M , . 
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wherein: 



F(n) is the nth value of the single sided scaled FFT spectrum; and 
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W represents a window function, wherein the window function is shifted 
by a value of the tone frequency variable fj. 

42. The method of claim 36, wherein said frequency transform comprises a 
windowed Fast Fourier Transform (FFT), and wherein said generating a frequency 
transform of the samples comprises generating a power spectrum of the samples; 

wherein the first expression and the second expression have the form: 

YJ\nn)W{f„ -/.)!) Sr(«)|xK(/„ 

_M _ k+\ 

Jfc+A/ Jt+M-l 

Y,\W{f„-L)\ Y.\Wifn-fi)\ 

wherein: 

F(n) is the nth value of the single sided scaled FFT spectrum; and 
W represents a window function, wherein the window function is shifted 
by a value of the tone frequency variable f^. 

43. The method of claim 36, further comprising: 
storing the determined tone frequency value in a memory. 

44. The method of claim 36, further comprising: 
outputting the determined tone frequency value. 

45. The method of claim 36, further comprising: 

computing one or more of the amplitude and phase of the tone using the 
determined tone frequency value. 
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46. The method of claim 36, wherein the pluraUty of expressions include: 

a first expression comprising a sum of one or more numerator terms divided by a 
sum of one or more denominator terms; and 

a second expression comprising a sum of at least two niunerator terms divided by 
a sum of at least two denominator terms. 

47. The method of claim 36, wherein the plurahty of expressions include: 

a first expression comprising a sum of two or more numerator terms divided by a 
sum of two or more denominator terms; and 

a second expression comprising a sum of three or more numerator terms divided 
by a sum of three or more denominator terms. 

48. The method of claim 1 , fiirther comprising: 

determining a plurality of potential local peaks based on the fi-equency response 
data for the plurality of tones, wherein each potential local peak corresponds to a tone; 

wherein said one or more signal parameter estimates are determined based on the 
plurality of potential local peaks; 

49. The method of claim 1 , 

wherein said generating modified frequency response data for the respective tone 
comprises generating modified fi"equency response data of the potential local peak of the 
respective tone. 

50. The method of claim 49, wherein said generating modified fi-equency 
response data of the potential local peak comprises: 

determining two or more sample values proximate to the potential local peak; and 
for each of the other tones. 
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calculating an effect value for each of the two or more sample values using 
the estimated signal parameters for the tone; and 

subtracting the respective effect values from each of the two or more 
sample values to generate a corresponding two or more modified sample values, wherein 
5 said two or more modified sample values comprise said modified fi*equency response 
data. 

51. The method of claim 50, wherein said two or more sample values are 
consecutive sample values. 

10 

52. The method of claim 50, 

wherein said frequency transform comprises a windowed Fast Fourier Transform 
(FFT) with window fimction W(f)^ wherein / denotes frequency; 

wherein said estimated signal parameters for the tone comprise an amplitude Au a 
15 frequency /y, and a phase (j>u and wherein said calculating an effect value for each of the 
two or more sample values using the estimated signal parameters for the tone comprises 
calculating: 

dF{m) = A.e^'^WUm + /i) + Ae'^'Wif^ - ) , 

20 

for each sample m, wherein m denotes an index of the sample, and^ denotes the 
frequency bin of the sample; and 

wherein subtracting the respective effect value from each of the two or more 
sample values comprises subtracting a respective value of dF(m) from each of the two or 
25 more sample values. 
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53. The method of claim 1, wherein, in each iteration of said generating 
modified fi^equency response data, (d) further comprises: 

applying a single tone estimation method to the modified fi'equency response data, 
thereby generating a refined estimate of the one or more signal parameters of the 
5 respective tone; and 

updating said one or more signal parameter estimates with the refined estimate of 
the one or more signal parameters of the respective tone; 

wherein the final one or more signal parameter estimates for each of the plurality 
of tones are determined based on the refined estimates of the one or more signal 
10 parameters of the respective tone. 

54. The method of claim 1, fiirther comprising: 
for each respective tone of the plurality of tones: 

applying a single tone estimation method to the modified fi-equency 
15 response data, thereby generating a refined estimate of the one or more signal parameters 
of the respective tone; and 

updating said one or more signal parameter estimates with the refined 
estimate of the one or more signal parameters of the respective tone; 

wherein said applying and said updating are performed for each of the plurality of 
20 tones in an iterative fashion; and 

wherein the final one or more signal parameter estimates for each of the plurality 
of tones are determined based on the refined estimates of the one or more signal 
parameters of the respective tone. 

25 55. The method of claim 54, wherein said single tone estimation method 

comprises: 

identifying a fi-equency location proximate to an amplitude peak in the fi-equency 
transform, wherein the amplitude peak corresponds to the tone; 
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selecting two or more frequency bins proximate to the frequency location in the 
frequency transform; and 

determining a tone frequency value that minimizes a difference between at least a 
first expression and a second expression; 
5 wherein the first expression comprises a sum of two or more numerator 

terms divided by a sum of two or more denominator terms, wherein the first expression 
includes a tone frequency variable; 

wherein each numerator term and each denominator term corresponds to 
one of the frequency bins; and 
10 wherein the second expression comprises a sum of one or more of the 

numerator terms divided by a sum of one or more of the denominator terms, wherein the 
second expression includes the tone frequency variable; and 

wherein the tone frequency value comprises a frequency of the tone. 

15 56. The method of claim 55, wherein said determining a tone frequency value 

comprises: 

computing a plurality of differences between the first expression and the second 
expression for different respective tone frequency values of the tone frequency variable; 
and 

20 selecting the tone frequency value that produces a smallest difference; 

wherein said computing a plurality of differences comprises performing a 
Newton-Rhapson root finding method. 

25 57. A memory medium comprising program instructions for determining one 

or more signal parameters for a plurality of tones in an input signal, wherein the program 
instructions are executable by one or more processors to implement: 



Atty. Dkt. No.: 5150-62101 



Page 60 



Meyertons, Hood, Kivlin, Kowert & Goetzel PC 



a) receiving samples of the input signal, wherein the input signal includes a 
plurality of tones; 

b) computing a frequency transform of the input signal using said samples, 
wherein the frequency transform comprises frequency response data for the plurality of 

5 tones; 

c) determining a plurality of potential local peaks based on the frequency 
response data, wherein each potential local peak corresponds to a tone; 

d) determining one or more signal parameter estimates characterizing each of 
the plurality of tones based on said plurality of potential local peaks; 

10 e) storing said one or more signal parameter estimates for each of the 

plurality of tones; 

f) for each respective tone characterized by said one or more signal 
parameter estimates: 

fl) removing effects of other tones of the plurality of tones from the 
15 frequency response data of the potential local peak of the respective tone using the stored 
one or more signal parameter estimates of each of said other tones, thereby generating 
modified frequency response data for the respective tone; 

f2) applying a single tone estimation method to the modified 
frequency response data of the respective tone, thereby generating a refined estimate of 
20 the one or more signal parameters of the respective tone; and 

G) updating said one or more signal parameter estimates of the 
respective tone with said refined estimate; 

wherein, after said removing, said appljong, and said updating for each tone, the 
one or more signal parameter estimates comprise refined estimates for each tone; and 
25 g) performing step f) one or more times to generate final one or more signal 

parameter estimates for each of the plurality of tones. 

58. The memory medium of claim 57, 
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wherein the one or more signal parameter estimates characterizing the pluraUty of 
tones comprise a knowledge base of said signal parameter estimates; 

wherein said storing one or more signal parameter estimates comprises: 

storing a copy of said knowledge base, wherein said copy comprises a 
5 prior knowledge base; and 

wherein said updating said one or more signal parameter estimates with the 
refined estimate of the one or more signal parameters of the respective tone comprises: 

updating the knowledge base with said refined estimate of the one or more 
signal parameters of the tone, wherein, after said removing, said applying, and said 
10 updating for each tone, the knowledge base is comprised of said refined estimates for 
each tone. 

59. The memory medium of claim 58, wherein said performing step f) one or 
more times comprises performing step f) in an iterative manner until one or more of 1) 
15 the number of iterations exceeds a threshold number of iterations, or 2) the difference 
between values of the one or more signal parameters of said knowledge base and values 
of the one or more signal parameters of said prior knowledge base is less than a threshold 
value. 

20 60. The memory medium of claim 57, wherein said one or more signal 

parameter estimates comprise estimates for one or more of fi-equency, amplitude, and 
phase. 

61. The memory medium of claim 60, wherein said estimate of firequency 
25 comprises the fi*equency value of the determined potential local peak. 

62. The memory medium of claim 60, wherein said estimate of amplitude is 

zero. 
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63. The memory medium of claim 60, wherein said estimate of ampUtude is 
the value of the potential local peak. 

5 64. The memory medium of claim 60, wherein said initial estimate of phase is 

zero. 

65. The memory medium of claim 57, wherein said frequency transform 
comprises a windowed Fast Fourier Transform (FFT). 

10 

66. The memory medium of claim 57, wherein said storing said one or more 
signal parameter estimates for each of the plurality of tones comprises: 

sorting said signal parameter estimates for the plurality of tones by amplitude of 
said potential local peaks; 
15 wherein said removing effects, said applying a single tone estimation, and said 

updating said one or more signal parameter estimates are performed for each of said 
pluirality of tones in order of decreasing ampUtude of said potential local peaks. 

67. The memory medium of claim 57, wherein said removing effects of other 
20 tones from the frequency response data of the potential local peak of the respective tone 

comprises: 

determining two or more sample values in a neighborhood of the potential local 
peak; and 

for each of the other tones, 
25 computing an effect value for each of the two or more sample values using 

the estimated signal parameters for the tone; and 

subtracting the respective effect value from each of the two or more 
sample values to generate a corresponding two or more modified sample values, wherein 
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said two or more modified sample values comprise said modified firequency response 
data; and 

wherein said modified fi-equency response data are usable in applying said single 
tone estimation method. 

5 

68. The memory medium of claim 67, wherein said two or more sample 
values are consecutive sample values. 

69. The memory medium of claim 67, 

10 wherein said fi-equency transform comprises a windowed Fast Fourier Transform 

(FFT) with window fimction W(f), wherein / denotes fi-equency; 

wherein said estimated signal parameters for the tone comprise an ampUtude ^4;, a 
frequency /y, and a phase and wherein said computing an effect value for each of the 
two or more sample values using the estimated signal parameters for the tone comprises 

15 computing: 

dF{m) = Ae^'^W{f^-^f,)-^Ae^'^W{f^ 

for each sample m, wherein m denotes an index of the sample, and^^ denotes the 
20 frequency bin of the sample; and 

wherein subtracting the respective effect value from each of the two or more 
sample values comprises subtracting a respective value of dF(m) from each of the two or 
more sample values. 

25 70. The memory medium of claim 57, wherein said single tone estimation 

method comprises: 
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identifying a frequency location proximate to an amplitude peak in the frequency 
transform, wherein the amplitude peak corresponds to the tone; 

selecting two or more frequency bins proximate to the frequency location in the 
frequency transform; 

determining a tone frequency value that minimizes a difference between at least a 
first expression and a second expression, 

wherein the first expression comprises a sum of two or more numerator 
terms divided by a sum of two or more denominator terms, wherein the first expression 
includes a tone frequency variable; 

wherein each numerator term and each denominator term corresponds to 
one of the frequency bins; and 

wherein the second expression comprises a sum of one or more of the 
numerator terms divided by a sum of one or more of the denominator terms, wherein the 
second expression includes a tone frequency variable which represents a correct tone 
frequency value of the tone; and 

wherein the tone frequency value comprises a frequency of the tone. 

7 1 . The memory medium of claim 70, 

wherein the tone frequency variable represents a correct tone frequency value of 
the tone; and 

wherein the determined tone frequency value represents an approximation of the 
correct tone frequency value. 

72. The memory medium of claim 70, 

wherein a ratio of each numerator term and its corresponding denominator term 
represent a complex ampUtude of the tone at a respective bin. 
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73. The memory medium of claim 70, wherein the first expression and the 
second expression have the form: 
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wherein: 

F(n) is the nth value of the single sided scaled FFT spectrum; and 
W represents a window function, wherein the window function is shifted 
by a value of the tone firequency variable fj. 

74. The memory medium of claim 70, 

wherein said generating a frequency transform of the samples comprises 
generating a power spectrum of the samples; 

wherein the first expression and the second expression have the form: 
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wherein: 

F(n) is the nth value of the single sided scaled FFT spectrum; and 
20 W represents a window function, wherein the window function is shifted 

by a value of the tone fi-equency variable fj. 

75. The memory medium of claim 57, further comprising: 
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storing the final signal parameter estimates for each of the plurality of tones. 



76. The memory medium of claim 57, further comprising: 

outputting the final signal parameter estimates for each of the plurality of tones. 

5 

77. The memory medium of claim 57, wherein the fi*equency transform 
comprises a fi-equency transform array for the plurality of tones, and wherein said single 
tone estimation method comprises: 

identifying two or more frequency bins proximate to a first magnitude peak in the 
1 0 frequency transform array; 

determining a tone frequency value that minimizes a difference between a first 
expression and a second expression; and 

selecting a tone frequency value that makes a plurality of expressions most nearly 

equal; 

15 wherein each of the plurality of expressions comprises a sum of one or 

more numerator terms divided by a sum of one or more denominator terms, wherein each 

of the plurality of expressions includes a tone frequency variable; 

wherein each numerator term and each denominator term corresponds to 

one of the frequency bins; and 
20 wherein a ratio of each numerator term and its corresponding denominator 

term represent a complex amplitude of the tone. 



78. A tone detection system, comprising: 
25 an input for receiving samples of an input signal, wherein the input signal 

includes a plurality of tones; 
a processor; and 
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a memory mediimi coupled to the processor, wherein the memory medium stores 
a tone detection software program for detecting the plurality of tones present in the input 
signal; 

wherein the processor is operable to execute the tone detection software program 
5 to implement: 

a) generating a frequency transform of the samples, wherein the frequency 
transform comprises frequency response data for the plurality of tones; 

b) determining a plurality of potential local peaks based on the frequency 
response data, wherein each potential local peak corresponds to a tone; 

10 c) determining one or more signal parameter estimates characterizing each of 

the plurality of tones based on said plurality of potential local peaks; 

d) storing said one or more signal parameter estimates for each of the 
plurality of tones; 

e) for each respective tone characterized by said one or more signal 
1 5 parameter estimates, 

el) removing effects of other tones of the plurality of tones from the 

frequency response data of the potential local peak of the respective tone using said one 

or more signal parameter estimates of each of said other tones to generate modified 

frequency response data for the respective tone; 
20 e2) applying a single tone estimation method to the modified 

frequency response data, thereby generating a refined estimate of the one or more signal 

parameters of the respective tone; and 

e3) updating said one or more signal parameter estimates with the 

refined estimate of the one or more signal parameters of the respective tone; 
25 wherein, after said removing, said applying, and said updating for each tone, the 

one or more signal parameter estimates comprise refined estimates for each tone; and 

f) performing step e) one or more times to generate final one or more signal 
parameter estimates for each of the plurality of tones. 
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79. The tone detection system of claim 78, wherein said removing effects of 
other tones from the frequency response data of the potential local peak of the respective 
tone comprises: 

5 determining two or more sample values proximate to the potential local peak; 

for each of the other tones, 

calculating an effect value for each of the two or more sample values using 
the estimated signal parameters for the tone; and 

subtracting the respective effect value from each of the two or more 
10 sample values to generate a corresponding two or more modified sample values, wherein 
said two or more modified sample values comprise said modified frequency response 
data; and 

wherein said modified frequency response data are usable in applying said single 
tone estimation method. 

15 

80. The tone detection system of claim 79, wherein said two or more sample 
values are consecutive sample values. 

8 1 . The tone detection system of claim 79, 

20 wherein said frequency transform comprises a windowed Fast Fourier Transform 

(FFT) with window ftmction W(j)^ wherein / denotes frequency; 

wherein said estimated signal parameters for the tone comprise an amplitude a 
frequency fu and a phase and wherein said calculating an effect value for each of the 
two or more sample values using the estimated signal parameters for the tone comprises 

25 calculating: 

dF{m) = A,e-^^W{f„ +f,) + - /,) , 
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for each sample w, wherein m denotes an index of the sample, and fm denotes the 
frequency bin of the sample; and 

wherein subtracting the respective effect value from each of the two or more 
5 sample values comprises subtracting a respective value of dF(m) from each of the two or 
more sample values. 

82, The tone detection system of claim 78, wherein said single tone estimation 
method comprises: 

10 identifying an amplitude peak in the frequency transform, wherein the amplitude 

peak corresponds to the tone; 

selecting two or more frequency bins in a neighborhood of the amplitude peak in 
the frequency transform; and 

determining a tone frequency value that minimizes a difference between at least a 
1 5 first expression and a second expression; 

wherein the first expression comprises a sum of two or more numerator 
terms divided by a sum of two or more denominator terms, wherein the first expression 
includes a tone frequency variable; 

wherein each numerator term and its corresponding denominator term 
20 corresponds to one of the frequency bins; and 

wherein the second expression comprises a sum of one or more of the 
numerator terms divided by a sum of one or more of the denominator terms, wherein the 
second expression includes the tone frequency variable; and 

wherein the first expression is approximately equivalent to the second 
25 expression when the correct tone frequency value is used for the tone frequency variable 
in the first and second expressions; and 

wherein the tone frequency value comprises a frequency of the tone. 
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83. The tone detection system of claim 82, 

wherein the tone frequency variable represents a correct tone frequency value of 
the tone; and 

wherein the determined tone frequency value represents an approximation of the 
5 correct tone frequency value. 

84. The tone detection system of claim 82, 

wherein a ratio of each numerator term and its corresponding denominator term 
represent an amplitude of the tone at a respective bin. 

10 

85 . The tone detection system of claim 82, 

wherein a ratio of each numerator term and its corresponding denominator term 
represent a complex amplitude of the tone at a respective bin. 



15 

86. A method for determining signal parameters for a plurality of tones in an 
input signal, the method comprising: 

a) receiving samples of the input signal, wherein the input signal includes a 
plurality of tones; 

20 b) applying a windowed FFT to the signal to generate a frequency transform 

of the signal, wherein the frequency transform comprises frequency response information; 

c) determining a plurality of potential local peaks based on the frequency 
response information, wherein each potential local peak corresponds to a tone; 

d) building a knowledge base of signal parameters characterizing said 
25 plurality of tones by determining initial estimates of signal parameters for each of the 

plurality of tones; 

e) storing a copy of said knowledge base of tone parameters, wherein said 
copy comprises a prior knowledge base; 
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f) for each respective tone characterized in said knowledge base, 

fl) removing effects of the others of the pluraUty of tones from said 
estimates of the signal parameters using said knowledge base to generate modified 
frequency response information for the respective tone; 
5 f2) applying a single tone estimation method to the modified 

frequency response information, thereby generating a refined estimate of the signal 
parameters of the tone; and 

f3) updating said knowledge base with said refined estimate of the 
signal parameters of the tone; 
10 wherein, after said removing, said applying, and said updating for each tone, the 

knowledge base is comprised of said refined estimates for each tone; and 

g) performing steps e) through i) in an iterative manner until 1) the number of 
iterations exceeds a threshold number of iterations, or 2) the difference between values of 
the signal parameters of said knowledge base and values of the signal parameters of said 

1 5 prior knowledge base is less than a threshold value. 



87. A memory medium comprising program instructions for determining one 
or more signal parameters for a plurality of tones in an input signal, wherein the program 
20 instructions are executable by one or more processors to implement: 

a) receiving the input signal, wherein the input signal includes a pluraUty of 

tones; 

b) computing a frequency transform of the input signal, wherein the 
frequency transform comprises frequency response data for the plurality of tones; 

25 c) determining one or more signal parameter estimates characterizing each of 

the plurality of tones based on said frequency response data; 

d) for each respective tone of the plurality of tones, generating modified 
frequency response data for the respective tone using the one or more signal parameter 



Atty. Dkt. No.: 5150-62101 



Page 72 Meyertons, Hood, Kivlin, Kowert & Goetzel PC 



estimates of one or more of the other tones, wherein said generating the modified 
firequency response data for the respective tone removes effects of other tones of the 
pluraUty of tones; 

wherein step d) is performed in an iterative fashion; and 
5 e) determining final one or more signal parameter estimates for each of the 

pliirality of tones based on the modified fi-equency response data of the plurality of tones. 
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